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Abstract 

We consider the scalar sector of the most general renormalizable two-Higgs-doublet 
model at non-zero temperature. We calculate the largest finite temperature corrections 
to the free-energy density and study thermal evolution of the ground state. Within the 
approximation chosen, we establish all possible sequences of thermal phase transitions 
and study their relation with the symmetries of the model. We show, in particular, 
that a charge-breaking or a CP- violating vacuum can arise at intermediate stages of 
thermal evolution, and that the first-order phase transition is associated with a discrete 
symmetry of the potential, but not of the entire scalar lagrangian. 

PACS numbers: 12.60.Fr; 14.80.Cp 

1 Introduction 

1.1 Why general 2HDM? 

The two-Higgs-doublet model (2HDM) is one of the simplest beyond-the-Standard-Model ex- 
tensions of the Higgs mechanism of the electroweak symmetry breaking, [TJ [2]. In this model 
one introduces two doublets of Higgs fields <p\ and 02, which interact with the matter fields 
and also self-interact via an appropriate Higgs potential. 

When constructing the Higgs potential of 2HDM, one encounters many free parameters, 
which are not constrained by experiment. Extensive studies conducted over past decades 
have shown that playing even with a small subset of these free parameters one can get a rich 
spectrum of models with different and quite remarkable phenomenological consequences (see 
examples in [2]). 

Still, there is something to worry about in this zoo of particular versions of 2HDM. First, 
it has been noted that in certain cases seemingly different 2HDM models lead to very similar, 
and sometimes identical, physics. These sporadic observations evolved later into the concept of 
reparametrization-invariance (or Higgs-basis-invariance) of the model. Understanding began 
to grow that it is the set of basis-invariant features of the theory, and not the entire set of 
free parameters, that really shapes the phenomenology. Papers [HI H] were driven precisely 
by the search for an efficient basis-invariant description of the general 2HDM. However, a 
method how to extract these basis-invariant features from the conventional parametrization 



of the potential remained elusive for a long time; search for basis-invariant conditions for the 
explicit CP- violation conducted in [5] nicely illustrates this kind of difficulties. 

Second, starting from the very first paper on this subject, [ili], there has always been a gen- 
eral understanding that extra symmetries of the 2HDM lagrangian have a profound effect on 
the resulting phenomenology. Indeed, majority of papers on 2HDM were centered around the 
concepts of CP- violation in the scalar sector, or additional Z2 or U(l) symmetries of the poten- 
tial. It has been noted that such symmetries can be implemented in seemingly different ways, 
still leading to similar physics, see e.g. recent analysis of generalized CP-transformations in 
[6] . However, a complete list of possible symmetries and their phenomenological consequences 
was very difficult to establish in a traditional approach and was missing until very recently. 

These issues left many researchers wondering if the most general treatment of 2HDM can 
be carried through at all. If that were possible, one could generate the exhaustive list of 
possibilities offered by two Higgs doublets, gain a deeper understanding of various particular 
versions of 2HDM and clarify relations among them. This would also become a useful play- 
ground for model-builders, who wanted to construct even more complicated Higgs sectors with 
predefined properties. 

In addition, one should not forget the Minimal Supersymmetric Standard Model (MSSM), 
which has always been an important motivation for the study of 2HDM. Although the Higgs 
potential in MSSM is very constrained by the supersymmetry, with several relations among 
its parameters, supersymmetry-breaking loop corrections will lift them. Thus, as stated in [3], 
the Higgs sector of the (broken) MSSM can be described by an effective field theory consisting 
of the most general 2HDM. 

Thus, the most general 2HDM with no a priori restriction on its free parameters should 
be viewed as a useful tool rather than a viable approximation to the reality, and it is definitely 
worth studying in as much detail as possible. 

1.2 Problems with the general 2HDM 

The problem with the general 2HDM is that it cannot be worked out with straightforward 
algebra. The obstacle arises at the very first step: when minimizing the Higgs potential, one 
arrives at algebraic equations of sixth order, which cannot be solved in the general case. In 
this situation, any method that would give any non-trivial insight into the most general 2HDM 
is very welcome. 

Recently, in [?l [9] a geometric approach was developed that allows one to analyze many 
features of the most general 2HDM circumventing the algebraic difficulty. This approach led to 
proof of various coexistence theorems, classification of the number and symmetry properties of 
the minima of the potential, and eventually to the complete description of the phase diagram 
of the scalar sector of 2HDM. All this was done without explicitly computing the exact position 
of the global minimum of the potential. 

Those results constituted the very first step towards the phenomenology of the most general 
2HDM. In this paper we make another step: we analyze how the ground state of the model 
changes at non-zero temperature. We calculate thermal corrections to the Higgs potential in 
the lowest non-trivial approximation (first order in Aj and enhanced by T 2 -factor), and study 
the thermal phase transitions that take place. It is remarkable that this, too, can be done 
without explicitly minimizing the free-energy density. 

We stress once again that it is not our aim here to provide a phenomenologically viable 



2 



approximations of reality. What we do here is we develop the only formalism available so 
far that works for the most general 2HDM, hoping that its results will be used later in more 
phenomenology-oriented studies. 

Thermal field theory is a very developed branch of QFT. Thermal phase transitions 
in multicomponent scalar field theories have been, of course, studied previously, see e.g. 
PH EH H21 [131 El H5]- A detailed analysis of the Higgs potential at finite temperature in 
2HDM was conducted for example in [HI [15]. We also draw reader's attention to the paper 
[16] where an attempt was made to describe the phase diagram of the 2HDM and study the 
phase transitions possible, which we find unsatisfactory and do not support some of its claims, 
see details below. However, the authors of all of these works, even though starting with rather 
general 2HDM, quickly focus on some rather simple cases, which could be analyzed by simple 
algebra. Thus, the most general 2HDM at finite temperature has never been analyzed up to 
now in any approximation. 

The structure of the paper is the following. In Section 2 we briefly review the geometric 
approach to the most general 2HDM. Then, in Section 3 we describe the phase diagram of the 
scalar part of the model. In Section 4 we analyze thermal evolution of the ground state and 
list possible sequences of phase transitions. In Section 5 we discuss how symmetries of 2HDM 
are related with possible thermal phase transitions. Discussion and conclusions are presented 
in Sections 6 and 7. 



2 Brief overview of the formalism 



In this work we focus on the scalar sector of 2HDM. The Higgs potential of the most general 
renormalizable 2HDM, Vh — V2 + V4, is conventionally parametrized as 
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It contains 14 free parameters: real mf 1( m\ 2) Ai, A2, A3, A4 and complex mf 2 , A5, A6, A7. As a 
result the 2HDM a very rich phenomenology even at tree level. However, different points in this 
14- dimensional parameter space do not necessarily correspond to distinct physics. Indeed, if we 
perform any linear transformation between doublets 0i and 02, we arrive at the same generic 
potential with redefined coefficients, which still corresponds to the same set of observables. 
Thus, the problem of minimizing the Higgs potential has a reparametrization invariance, with 
reparametrization group GL(2,C). This property is an extension of an earlier identification 
of rephasing U(l), and unitary reparametrization U(2), invariance of the model, [31 141 [T§1 [2D]; 
it was developed in detail in [U] and helped prove several properties of the general 2HDM. 
In this Section we will briefly review this formalism. 

It starts with the observation of the internal Minkowski-space structure behind 2HDM. 
We first introduce the four-vector r M = (ro, r«) = ($^<3>, where $ = 02) T is a 

2-dimensional vector of Higgs doublets and a % are the Pauli matrices. This four vector is 
gauge invariant and parametrizes the gauge orbits in the space of the Higgs fields. 
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The general reparametrization group GL(2, C) can be written as C* x SL(2, C), where C* 
is the group of simultaneous multiplication of both fa with the same complex number, while 
SL(2, C) is the special linear transformation group containing all unit-determinant transfor- 
mation matrices. Multiplication of both doublets by the same number gives a freedom to 
arbitrarily rescale r M , however it does not modify the structure of the Higgs potential. The 
special linear group, on the contrary, leads to non-trivial transformations of the Higgs poten- 
tial, and in subsequent analysis we focus on it. 

Transformations of $ under the special linear group SL(2, C) corresponds to the SO(l, 3) 
transformations of r M , equipping the gauge orbit space with the Minkowski-space structure. 
It follows from the definition of r M that 



r = + (4>\fa) > 



r^ry, = 4 



(<t>\fa)(4fa) - (cf>\fa) (4fa) 



> 0. 



Note that r M r M > is possible thanks to the fact that fa and fa are themselves electroweak 
doublets, and the last line of <^ is just the Schwarz inequality written for fa. Thus, the 
physically realizable vectors r M populate not the entire 1 + 3-dimensional Minkowski space, 
but the future lightcone. The Higgs potential in the r M -space can be written in a very compact 
form: 

K = -M/ + ^VV. (3) 

Here the four- vector is built from parameters m?- in ([I]), while the symmetric four-tensor 
is constructed from the quartic coefficients Aj. Their explicit expressions as well as some 
properties can be found in [HI E]- Here we just note the most important property of A^: it 
can always be diagonalized by a specific 50(1,3) transformation of the r^-space, and after 
diagonalization it takes form 

A^ = diag(A , -A 1; -A 2 , -A 3 ) with A > 0, A > A; , i = 1, 2, 3 , (4) 

where the inequalities among the eigenvalues result from the positivity constraint on the 
potential. The minus signs in front of the "space-like" eigenvalues arise from the pseudo- 
euclidean metric in the orbits space. 

It is known that the potential ([3]) can have three types of minima: (i) the electroweak 
(EW) conserving, (ii) the EW-breaking but charge conserving (i.e. neutral), and (iii) the EW- 
and charge-breaking ones. One can use the v.e.v.s of the two doublets (fa) to construct (r M ). 
Then, the three type of minima correspond to: (i) (r^) = (the apex of the forward lightcone 
LC+), (ii) (r^) ^ but (r^)(r M ) = (the surface of LC+), (iii) (r^) ^ and (r")<r M ) > (the 
interior of LC + ). 

The position of the charge-breaking extremum (r u ) c h is given by the following coupled 
equations: 

A^(r v ) ch = M», (5) 

If A^ is not singular, a solution of this system always exists and is unique. However solutions 
realizable in the 2HDM, i.e. such that (r u ) c h lies inside the forward lightcone, exist only if 
M M lies inside some cone with the apex at the origin. If this is the case, then an additional 
requirement (Aj < 0, i — 1, 2, 3, i.e. the tensor A^ IU is positive-definite in the entire space 
of non-zero vectors r M ) guarantees that the charge-breaking extremum is a minimum, [8]. If, 
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on the contrary, at least one of A, is positive, then the charge-breaking extremum is a saddle 
point. 

The positions of all neutral extrema (r^) are the solutions of the following coupled equa- 
tions: 

A""<r„) - C<0 = M" , (r")(r M ) = , (6) 

where ( is a Lagrange multiplier. This system can have up to six solutions. To distunguish 
a (local) minimum from other types of extrema one needs to use more involved tools, for 
example, geometric methods developed in [H [9]. These methods allowed to prove that the 
generic potential has at most two local minima in the entire orbit space (not just on the surface 
of the cone), while the other extrema are saddle points. 

Finally, the Higgs potential (CD) must be accompanied by the kinetic term, which can also 
be written covariantly: 

K = K lif f, p^=(D a ^a»(D a $), (7) 

where D a is the covariant derivative, a denotes the usual space-time coordinates, while /i, as 
before, denotes the coordinate in the orbit space. Note that reparametrization transformation 
properties of p M are the same as r M . In the default frame, = (1, 0, 0, 0). Upon an 
5*0(1, 3) transformation, acquires non-zero "space-like" coordinates, however the condition 
K^K^ = 1 is always satisfied (if we considered GL(2, C) transformations of the fields instead 
of SL(2,C), this would not hold anymore; however we do not use the rescaling subgroup of 
GL(2, C) just because the SL(2, C) group itself is sufficient in order to describe the phases 
and the sequences of possible phase transitions of 2HDM). The four- vector is not involved 
in the search for the minimum of the potential, however it affects the mass matrix at this 
minimum. 

This generalized kinetic term effectively incorporates the non-diagonal kinetic term, which, 
as was argued in [21] , must be introduced in the initial lagrangian to restore renormalizability 
of the model. 

3 The phase diagram of the scalar sector of 2HDM 

Before studying thermal evolution of the 2HDM vacuum, let us explicitly describe the phase 
diagram of the zero-temperature 2HDM. Here we use results from [9[ [22] but slightly refor- 
mulate them in anticipation of the study of thermal evolution. 

Note that whenever we say "phase transition" in this Section, we imply zero-temperature 
phase transition, that is, transitions from one phase to another upon continous change of the 
parameters of the potential. We will classify such zero-temperature phase transitions as the 
"second-order" and "first-order" ones according to whether the order parameter (the position 
of the global minimum of the potential) changes continuously or not. 

We first switch to the A^- diagonal frame. We will focus here on the generic situation, in 
which all the "space-like" eigenvalues Aj are different. The phases can be then represented as 
regions in the M M -space in this frame. 

The structure of the phase diagram depends on the sign of M . 

3.1 The lower half-space 

In the lower half-space of M M , i.e. with Mq < 0, the phase diagram is extremely simple. 
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• If M^M^ > 0, the only extremum of the potential is the minimum at the origin. Thus, 
the past lightcone in the M^-space is the EW-symmetric phase. 

• If M^M^ < 0, i.e. M M lies outside the past lightcone but still in the lower half-space, 
then there exists a unique non-zero neutral extremum, which is necessarily the global 
minimum. 

Thus, in the lower half-plane we have only one surface of phase transitions, the past lightcone, 
at which the EW-breaking/restoring phase transition takes place. 

3.2 The upper half-space 

For M > 0, the phase diagram is much richer. To describe it, let us introduce the 3- vector 
trij, i = 1,2, 3: 

= t^(M 1; M 2 , M 3 ) (8) 

and show the phase diagram in the mj-space. 

There are two generic cases to consider, which are shown in Fig. [TJ 




Figure 1: The phase diagram in the trij-space. Left: all Aj < 0; the ellipsoid separates the 
charge-breaking and the neutral vacua. Right A 2 is positive and greater than A x and A 3 . 
Inside the dark ellipse the discrete symmetry of the potential V is spontaneously broken. 



3.2.1 All A, < 

The condition Ai, A 2 , A3 < is the necessary condition for the charge- violating minimum to 
exist. Whether this minimum is indeed realized, depends on m,, i.e. on the point of the phase 
diagram. Indeed, in the A M „- diagonal frame, Eq. ([5]) can be rewritten as 

A (r ) = M , A i (r i )=M i . (9) 

Since the solution (r^) must lie inside the forward lightcone, one obtains that the charge- 
breaking phase in the phase diagrams lies inside the ellipsoid 

r 4 + V 4 + ^<l, where a, = JM , i = 1, 2, 3 , (10) 
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shown in Fig. [TJ left. 

If ttlj lies outside ellipsoid (jlOp . then the vacuum is neutral. The surface of the ellipsoid is 
the locus of the critical points of the phase diagram, at which the second order charge-breaking 
or charge-restoring phase transition takes place. 



3.2.2 At least one A; > 

If at least one Aj is positive, then the minimum always corresponds to a neutral vacuum and 
its position satisfies equations Writing (r M ) as r (l, ni, n 2 , n 3 ) with |n| = 1, one can 
eliminate the Lagrange multiplier from ([6]) to arrive at the following system of equations: 

[M -(A -A i )r Q ]n i = M i , i = 1,2,3. (11) 

The requirement that n is a unit vector can be written as 

^-^T2 = 1 - < 12 > 

• (l - ^) 

This is a sixth order equation for ro- 

In a general case, the system cannot be solved exactly. However, the geometric approach 
of [H [9] still allows one to understand the structure of the phase diagram. 

First, if irij lies on one of the principal planes, then the potential has an additional discrete 
symmetry. The solutions of (FIT]) can either conserve or violate this symmetry. The necessary 
and sufficient condition for violation of this symmetry, e.g. in the case m.2 = 0, are 

m? ml A 2 — Aj , . 

Thus, symmetry-violating extrema appear, if rrij lies inside a specific ellipse on one of the 
principal planes. 

In [SI [22] it was proved that these symmetry- violating extrema are minima, if and only if 
the corresponding eigenvalue of (A 2 for m 2 = 0) is positive and is the biggest of all Aj. So, 
when constructing the phase diagram, we should first identify the largest eigenvalue and then 
consider only the ellipse that lies in the plane orthogonal to the corresponding eigenvector, as 
it is shown in Fig. [TJ right. For points lying strictly inside the ellipse two different degenerate 
minima exist. In both of them the discrete symmetry is spontaneously broker^]. For points 
in the (mi, 1113) plane outside the ellipse only one minimum exists and the discrete symmetry 
is preserved. So, if rrij lies strictly on the plane and moves from outside into the ellipse, a 
symmetry-breaking second-order phase transition takes place. 

The points just above (rri2 > 0) or just below (rri2 < 0) the ellipse have two minima at 
different depths. So, if rrij lies above the ellipse and moves through it downwards, the relative 
depth between the two minima changes sign, and a first-order phase transition takes placed). 
Thus, the interior of the ellipse is the locus of the first-order phase transitions, while its 
boundary is the locus of the second-order phase transitions. 



With an appropriately symmetric kinetic term K^, this extends to the symmetry of the whole scalar 
lagrangian. The case depicted in Fig. [TJ right, corresponds then to spontaneous CP- violation. 

2 We consider here only the phase diagram at equilibrium and do not discuss the actual kinetics of this 
phase transition. 
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Let us also discuss the phase diagram of 2HDM obtained in [16]. The authors of that paper 
used not the most general, but a rather restricted version of the 2HDM potential, with only 8 
non-zero free parameters out of 14. The parameters \$, Xj and m\ 2 were set to zero, so that 
the potential from the start was (Z2) 2 -symmetric. This choice authomatically excludes any 
CP- violation, which is the origin of the authors' surprising assertion that the relative phase 
between the two v.e.v.'s "is not demonstrated in the particle properties". 

The potential chosen in [16] was sufficiently simple to allow for straightforward minimiza- 
tion. Six phases were identified and classified according to the residual symmetry and the 
"number of condensates", i.e. whether ((j)\<f>i) and {4>\<t>2) ar e non-zero. We find this picture 
misleading for the following reason. If one focuses on the scalar part of the theory only, then all 
the observables, including classification of the phases according to the properties of the ground 
state, must remain invariant under the group of the most general reparametrization transfor- 
mations of the scalar lagrangian. However, the classification of [16] is not reparametrization 
invariant. Phases II, III and IV are all related by reparametrization transformations and cor- 
respond just to the neutral vacuum, while phases V and VI correspond to the charge-breaking 
vacuum. This is not surprising since the form of the 2HDM potential used in [16] was rather 
restrictive and obscured these relations. Therefore, some of their phases do indeed correspond 
to truly different 2HDM phases: the governing principle is not the number of condensates, but 
the symmetries of the lagrangian and of the ground state. 



4 Thermal evolution of the ground state 
4.1 Temperature dependence of the mass terms 

Let us now consider the thermal evolution of the ground state of the model. As said in the 
Introduction, our motivation is not to provide a phenomenologically accurate description of 
thermal 2HDM, but rather to extend the developed formalism, which effectively works for 
the most general 2HDM, to non-zero temperatures. To this end, we stick to the following 
approximations: 

• we consider only the scalar sector of 2HDM and neglect temperature corrections from 
gauge bosons and the fermions, 

• we consider only the largest one- loop high-temperature corrections oc A;T 2 , 

• although the results can be trusted only at temperatures much larger than the mass 
parameters, we still use them to discuss the lower temperature region, T ~ m. 

Although in this region our results must be strongly corrected by terms which we neglect, we 
still think our discussion provided a first useful glimpse at the typical patterns of phase tran- 
sitions one can expect in 2HDM. It is possible that the approach used here can be developed 
further to remove some of these assumptions. 

The largest thermal correction to the scalar 4 -theory comes from the hard thermal one- 
loop self-energy diagrams, see e.g. [23]. In the one-component real scalar field theory with 
negligible zero-temperature mass and interaction term A</> 4 /4!, the thermal correction is 

5m 2 = . 
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In 2HDM, we have similar diagrams with various fields in the loop. In order to extract the 
largest, oc A«T 2 correction, we neglect mass terms in the propagators in loops, so that the 
scalar loops generate only a universal factor. Thus, the temperature corrections to m 2 - differ 
only in vertex factors, and the straightforward calculation gives: 

2 irr\ 2 / n \ rp2 3Ai + 2A 3 + A 4 



m 2 u (T) = miM-T 
m 2 22 (T) = m 2 22 (0)-T 



2 /-;m _ 2 ,ru -i--2 3A 2 + 2A 3 + A 4 



6 

ml 2 (T) = m 2 2 (0)-T 2 -^±^. (14) 



As we already said, strictly speaking, these equations are reliable only at temperatures much 
larger then masses of the particles in loops. At lower temperatures, the mass parameters 
cannot be neglected anymore, which significantly complicates the problem. However, the 
present study is only the first approximation to a more complete analysis of thermal 2HDM, 
so we adopt (fT4l) at all temperatures and follow its physical consequences. 

In the Minkowski-space formalism, one can rewrite (|14l) as the temperature dependence of 
the SO(l, 3)-covariant four-vector M M . Indeed, in the Appendix we show that 

M"(T) = AP(0) -c^-T 2 , c, = -^ (6A^ - <^TrA) K v , (15) 

where TrA = A^g^ = Aq + Ai + A 2 + A3 = A3 — A4. The four- vector of kinetic terms K v 
appears in this result due to the presence of the loop propagator. Note that K v does not 
receive any temperature corrections within this approximation. 

Expression (TT5T) can be used in any frame. In particular, in the original frame it reproduces 
(114)) . while in the A^„- diagonal frame it gives: 

6A - TrA 5A - Ai - A 2 - A 3 ^ ^ 
c o = ^ K = K > , (16) 

6Ai - TrA . , 

c i = —L— Ki t = l,2,3. (17) 

Note that can have any sign depending on the values of Aj. 

One can ask whether there is any constraint on the possible values of q. The answer is in 
the affirmative. Since K^K^ = 1, one obtains from (1151) : 

v c? / 6Ap - TrA V 2 _ v Kf 

\<c>\eA^A) ~\lq <l - (18) 

Thus, possible Cj/c lie inside the ellipsoid, whose semiaxes are always smaller than 5. 

Consider next M (T) in this frame: M (T) = M — CqT 2 . If M Q < 0, it further decreases 
with the increasing temperature without changing its sign. Thus, the only phase transition 
one can have in this case is the EW-restoring phase transition. 

If the initial Mq > 0, it will change its sign at 



T = T* = \I^. (19) 
c 
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Then one can separate thermal evolution of the ground state into two stages, which can be 
dubbed the "low-temperature" (T < T*) and the "high-temperature" (T > T*) stages. Again, 
it is only the low-temperature stage that non-trivial phase transitions can take place, while 
during the high-temperature stage one can observe only the EW-restoring phase transition. 
Focusing on the low-temperature stage, let us define m»(T) similarly to (jHJ): 

Mi{T) (M x -c x T 2 M 2 -c 2 T 2 M a 



tm(T) 



If we introduce x 



M (T) VMo - c T2' M - c T2' M - c T^ ' 
T 2 /(T 2 — T 2 ), so that x = at zero temperature, then 



(20) 



m,- 



1+1- 



c mi 



X 



m,(T) 



if m, ; 7^ . 



if m, = . 



Co 



(21) 



Thus, the temperature evolution rrij(T) is represented on the phase diagram by a straight ray 
starting from ttij at T = and going to infinity as T — > T*. 

It is easy to see that this result is not specific for the quadratic temperature dependence. 
It is based solely on the fact that the functional form of the temperature dependent correction 



to m. 2 ,- is the same. 



4.2 Sequences of thermal phase transitions 

The fact that thermal evolution of rrij is represented by a straight ray allows us to immediately 
classify possible sequences of phase transitions without any further calculation. 

Depending on the eigenvalues Aj, the sign of each of the components ttlj(T) in ( 1211) can 
either remain unchanged or can change only once at some temperature. 

If the signs of all rrij(T) do not change, the only phase transition one can have at growing 
temperature is the charge-restoring second-order phase transition, which takes place for nega- 
tive Aj and trij inside the ellipsoid, or the discrete symmetry restoring phase transition, when 
at least one Aj is positive and tTij lies inside the corresponding planar ellipse. 

If one or more rrij(T) change their signs, then one can have a richer spectrum of possibilities. 
On its way to infinity, the ray can go through various intermediate phases. In Fig. [2] we show 
examples of such paths for the two phase diagrams shown previously in Fig. [lj In Fig. [21 left, 
we show that the charge-breaking phase can be an intermediate phase of thermal evolution 
of the model. Fig. [2} middle, shows a path that leads to a first-order phase transition, while 
Fig. [2, right, demonstrates a planar ray that goes exactly through the ellipse and leads to a 
sequence of two second-order phase transitions: the symmetry-breaking and the symmetry- 
restoring ones. 

By combining various initial positions and directions of the thermal evolutions, one can 
easily generate a list of various sequences of phase transitions. However, from general geometric 
considerations (a straight line can intersect an ellipse or ellipsoid no more than twice), it follows 
that the longest chain of phase transitions consists of two phase transitions (followed by the 
electroweak phase transition at very high temperatures). Examples are shown in Fig. [21 

It can also happen that the ray barely touches but does not intersect the ellipsoid of Fig. [fl 
left, or the ellipse of Fig. [[J right. In such cases, the temperature evolution drives the system 
to a critical point and then returns back to the initial phase. 
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rri2 ITI2 




Figure 2: Examples of the temperature evolution trajectories (blue rays) for the two phase 
diagrams shown in Fig. [TJ In the rightmost figure the ray lies in the plane of the ellipse. 

5 Symmetries 

The scalar sector of the two-Higgs-doublet model can have various additional explicit symme- 
tries pi [19]. Suppose that 2HDM possesses some explicit symmetry at zero temperature. One 
can then ask: 

• does thermal evolution preserve this symmetry? 

• how are the possible sequences of phase transitions restricted by such a symmetry? 

A detailed analysis of possible symmetries of 2HDM was carried out in [9 J . It was stressed 
there that one should distinguish between symmetries of the potential and of the entire la- 
grangian, since they play different roles in the theory. A symmetry of the potential governs 
the vacuum structure and possible degeneracies, while a symmetry of the full Higgs lagrangian 
is needed to make the mass matrix symmetric, too. 

Symmetries are best explored in the A^-diagonal frame. If all A, are different, then in 
this frame the quartic part of the potential has symmetry group (Z 2 ) 3 generated by reflections 
of the three directions. If some of the spacelike components of JVf are zero in this frame, 
then the potential has some non-trivial explicit symmetry. If, in addition, the corresponding 
components of are also zero, then this symmetry extends to the entire lagrangian. 

Consider now Eq. ffl5|) . which gives M M (T). It is easy to see that if at T = there is 
a symmetry of the entire lagrangian (i.e. of A^, and K 11 ), then it will be preserved 
during the entire thermal evolution. If, instead, only the Higgs potential, but not the entire 
lagrangian, has some symmetry, then generically this symmetry is explicitly violated at non- 
zero temperatur^f). There is, however, an exceptional situation when Aj satisfy the relation 
TrA = 6Aj, so that the corresponding correction coefficient q = 0. In this case, the symmetry 
of the potential is preserved during the temperature evolution. 

Let us now take a closer look at the first-order phase transition, which happens when the 
evolution path crosses the ellipse of Fig. [2J middle. Exactly at the phase transition there are 
two separate degenerate minima of the Higgs potential. From Proposition 4 of [9] it follows 
that at this point (and only at this point!) a discrete symmetry of the potential emerges. At 

3 One must of course keep in mind that a symmetry of the potential not shared by the kinetic term is not 
stable in the perturbative series. A discussion involving such a symmetry is thus restricted to the tree-level 
approximation of the potential. 
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higher temperature this symmetry is lost again, therefore, it cannot be the symmetry of the 
entire lagrangian. Thus, we arrive at the following conclusion: 

The first-order phase transitions are associated with points of a discrete symmetry of the po- 
tential, but not of the entire lagrangian. 

One concludes that the symmetries which are often studied in literature, like CP ox the 
so-called ^-symmetry (which is in fact a (Z 2 ) 2 -symmetry, see [9]) cannot be associated with 
first-order phase transitions in 2HDM. An example of symmetry that can be associated with 
the first-order phase transition will be given in |17j . 

Of course, this result is derived only within the approximation chosen. Higher-order cor- 
rections to the free energy potential can turn some of the second-order phase transitions into 
the first-order ones. 

Let us also discuss the phase transitions identified in [16J. Similarly to our analysis, the 
authors of that paper assumed that the thermal evolution can still be described by the same 
fourth-order free energy density, but with the coefficients depending on temperature. They, 
however, did not calculate this dependence, thus overlooking a possibility of electroweak sym- 
metry non-restoration at large temperature. Among the other results, they find the possibility 
of first-order phase transition between the charge-breaking and the neutral vacua (in their no- 
tation, between phases IV and V). Our analysis does not leave any room for such a phase 
transition. The first-order phase transitions can arise only when two separate local minima 
coexist. However, we proved in [8] that disjoint charge-breaking and neutral minima never 
coexist in any 2HDM potential. Therefore, charge-breaking and charge-restoring transitions 
occur via a continuous change of the vacuum expectation values. The depth of the global 
minimum then changes smoothly, and there is no specific heat associated with such a phase 
transition. 

In fact, within the approximation chosen a first order phase transition can only take place 
between the same symmetry phases (isostructural phase transition). 

6 Discussion 

6.1 Non-restoration of the EW symmetry 

Since the seminal paper [TT] of Weinberg it is known that a symmetry, global or local, can 
remain broken even at very high temperatures. The simplest example of such a case is given 
by a 4 -model with two real scalars, [12] . This possibility is of course present in 2HDM, as 
can be seen directly from eq. ( fT4l) giving the temperature evolution of the mass terms m 2 ^ (T) . 
This happens, for example, when A 6 = A 7 = 0, and 3Ai + 2A 3 + A 4 and 3A 2 + 2A 3 + A 4 are of 
opposite signs. Together with the positivity conditions this typically implies that Ai A 2 (or 
vice versa), with 2A 3 + A4 sufficiently negative. 

In the Minkowskian formalism, this situation takes place, when M M (T) goes down in the 
phase diagram, but not sufficiently steeply to cross the past lightcone, i.e. if the four-vector 
c M defined in ( TTBT) is spacelike, c^c^ < 0: 

36(A 2 ) ^K^JC - 12TrA • A flI/ K fl K v + (TrA) 2 < . (22) 
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A necessary, but not sufficient, condition for this to take place is that there exists at least one 
negative Aj. 

Note also that the electroweak symmetry is the only symmetry that might remain broken 
at very large temperatures. All other symmetries are unaviodably restored, as it is obvious 
from the analysis of the phase diagram. This generalizes the result obtained in [13] that the 
spontaneously broken CP-symmetry is always restored at high temperatures in the 2HDM. 

6.2 Beyond the approximation chosen 

Within the approximation chosen, the above discussion gives an exhaustive list of sequences 
of phase transitions possible in 2HDM. However, various additional corrections will definitely 
modify this picture. 

Firstly, for temperatures T comparable with the mass parameters the loop integrals deter- 
mining the corrections to the mass terms of the potential are not universal. As a result, the 
functional T-dependence of these terms is more complicated and the thermal evolution will not 
follow the straight rays on the phase diagram. Second, there will be temperature corrections 
to Aj, which will result in temperature evolution of the phase diagram itself. Third, the other 
fields present in 2HDM will contribute to the effective potential and will affect evolution of 
the ground state. 

It remains to be seen whether the formalism used here can be extended further to cope 
with these terms, still preserving the generality of the model. 

7 Conclusions 

The two-Higgs-doublet model has a rich phase diagram even in the tree-level approximation. 
One can expect a similarly rich spectrum of phase transitions of 2HDM at non-zero tem- 
peratures. Unfortunately, analysis of thermal evolution of the most general 2HDM has been 
impossible for a long time due to algebraic complexity of the problem. In this paper, armed 
with the recently developed geometric approach to the most general renormalizable 2HDM, 
we get an insight into various properties of the model at non-zero temperature. 

We have computed the largest temperature corrections to the mass terms in a basis- 
invariant way and demostrated that, in the approximation adopted, the thermal evolution 
of the ground state is described by a straight ray on the fixed phase diagram. This allowed 
to establish all possible sequences of phase transitions within this approximation. Among 
the results that we have obtained are: the possibility to have spontaneously CP-violating 
or charge-violating ground states as intermediate stages of the thermal evolution, and the 
possibility to have first order phase transitions associated with a momentary restoration of a 
discrete symmetry of the potential, but not of the entire scalar lagrangian. 

These findings might have intriguing cosmological consequences, which will be studied in 
a more phenomenology-oriented paper [TT] . see also some preliminary discussion in [18j. The 
purpose of this work was just to show how to use the geometric approach to 2HDM to get 
insight into its thermal dynamics. 
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and useful comments, as well as to the referee for numerious suggestions for improvement of 
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A Derivation of the covariant temperature corrections 

Here we derive formula ffl5|) for the temperature corrections to the mass terms in the covariant 
form. 

We start with the Higgs lagrangian in the tensorial form: 

L = K a b (d^* a ){d^ b ) + - ZY^M^) . (23) 

Note also that to make the expression more symmetric, we assumed in (1231) that all terms are 
written separately, i.e. instead of As(0|0i)(0202) we write 

• • • + ^A 3 (0l0i)(A> 2 ) + \\M<p2)(4<Pi) + ■■■ (24) 

In contrast to the standard tensorial form, here indices a. . .d are double indices: 

a = {a, i} , a = 1,2, i =|, [ . 

The subindex a labels the doublets and i labels components of the doublet a. The values of 
Y a b and Z a b c d are insensitive to i: 

K a b = K {<*,i) {j3 . } = K « p . gi ; y\ = Y^\ m = Y a p ■ 5) , 

6 d = Z l f {l3J} {Sm} = Z p y s - d j d m . (25) 

Then, the 1-loop correction to the self-energy is given by the contraction of the quartic coupling 
with the propagator, which brings in K^ 1 : 

8Y\ = -C ■ {Z\ c d + Z a d \ + Z\ a d + Z\\) {K-\ d = -2C (Z\ c d + Z\\) (K' 1 )/ , (26) 

where C = T 2 /12. Note that summation over all possible ways of contracting indices au- 
tomatically takes into account configurations with identical fields. One can then sum over 
intra-doublet indices i, which gives: 

SY a p = -2C (2Z a ^ s + ZTfp) . (27) 

The extra factor 2 in the first term inside the brackets comes from 8j8j = 2. This expression 
is well-known, see e.g. Eq. (11) in [13]. 

The next step is to translate (I27p to the covariant notation. If A is a Hermitian 2x2 
matrix, then one can write 

A = aocr — a % o % = a^a^ , 
where a 11 = (a , a 1 ) anda^ = (er°, a 1 ). The inverse matrix can be written using a^ 1 = (o , —a 1 ): 

n 2 ' f* 
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If follows from definitions that Z a p' y s can be represented as 

Therefore, contraction of the first term inside brackets in (I2T1) gives 

2 • • ^ • Tr(o*aP) = 2A^K V , 

where for clarity we suppressed the remaining indices a and (3, and we also used K 2 = 1. To 
calculate the contraction of the second terms, we use identity 

l -{o p a p a v + a v a p a") = g pp a v + g vp a p - g pv a p , 

which gives 

hv^ ■ a p a p a u ■ K p = i(2A^ - TrAg^IC , 
with TrA = A liV g ,w = A3 — A4. Together, they make 

SY% = -C(6A^ - TrA^iTK)^ • 
Since tfF^ = 5M p {a^) a ^ we finally get 

5Mp = ~Tl (6A ^ ~~ ^ TrA) ^ ' (28) 
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